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We consider the evolution of the unstable periodic orbit structure of coupled chaotic systems. This 
involves the creation of a complicated set outside of the synchronization manifold (the emergent set). 
We quantitatively identify a critical transition point in its development (the decoherence transition) . 
For asymmetric systems we also describe a migration of unstable periodic orbits that is of central 
0^ ' importance in understanding these systems. Our framework provides an experimentally measurable 

0^ ' transition, even in situations where previously described bifurcation structures are inapplicable. 
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The idea that several subsystems, when interacting nonlinearly, collectively give rise to novel dynamics that are not 
obviously attributable to the individual component parts has been termed emergence @] . In this Letter we investigate 
such novel dynamics in systems of coupled chaotic maps, with an emphasis on systems of dissimilar components. When 
synchronized, the time evolution occurs on a restricted manifold (called the synchronization manifold) embedded in 
T—i , the full state space. As the degree of coupling is decreased to zero, the system gradually evolves into a completely 
£> ■ unsynchronized state in which all the degrees of freedom of the individual component maps are realized. At each 
extreme, the dynamics can be understood in terms of the components. In between, however, the situation is more 
complicated. 

Various transitions in this desynchronization process have been described in the literature p[-pT|. Much of this 
' earlier work depends on an invariant manifold that persists under decreased coupling due to system symmetries. Our 
0^ , approach is novel in that it does not refer to invariant manifolds or symmetries. Instead, we analyze the evolution of 
the periodic orbit structure as the coupling is decreased and synchronization breaks down. Our formalism is therefore 
applicable to a much larger class of coupled systems, in particular, those consisting of dissimilar components. We 
£>y report two main results. We describe the creation and evolution of a complicated set that develops outside of 
the synchronization manifold (the emergent set), and we quantitatively identify a critical transition point in its 
development (the decoherence transition). For asymmetric systems we also describe a migration of unstable periodic 
orbits that is of central importance in understanding these systems. Our framework is advantageous because it 
■ provides an experimentally measurable transition in situations where previously described bifurcation structures are 
?^ | inapplicable. 

Previous work has focused on the invariant dynamics in the synchronization manifold A4, which can easily be 
identified in coupled systems with symmetry (such as when two identical sub-systems are coupled together). On A4, 
the components evolve identically, and are said to exhibit identical synchrony [^). As the coupling decreases from 
a fully synchronized state, a bubbling bifurcation || occurs when an orbit within M (usually of low period | fl2|| ) 
loses transverse stability. In the presence of noise or small asymmetries, a typical trajectory quickly approaches and 
spends a long time in the vicinity of M., but makes occasional excursions. (If an appropriate attractor exists outside 
the synchronization manifold, this transition leads to the creation of riddled basins Q.) As the coupling is further 
decreased, the blowout bifurcation ||] is observed when M. itself becomes transversely unstable (on average). More 
recent work has described bifurcations that lead to the creation of periodic orbits off the synchronization manifold 
||; these may lead to the creation of chaotic attractors external to M. 0. Also, imperfect phase synchrony has been 
analyzed recently in terms of unstable periodic orbits Q , and synchrony transitions have been investigated in coupled 
lattices of identical maps Q . 

The concept of (differentiable) generalized synchrony (GS) |HJ[nj extends these ideas. GS relaxes the condition that 
the state variables evolve identically, and only requires that they be functionally related. As the coupling is reduced 
and GS breaks down, however, this function may become extremely complicated, and the identification of bubbling- 
type or blowout-type bifurcations is especially problematic. Thus, a more general description of the desynchronization 
process is needed. 

In the present work we extend the above description by considering the evolution of the unstable periodic orbits 
(UPOs) as the coupling is varied. We use the following two-dimensional, unidirectionally coupled system 
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f x -> /(a:) 

\ y -> c/(a;) + (1 - c)fli(y), 



(1) 



where / and g are chaotic maps and c is a scalar that describes the coupling. We emphasize that / and g need not 
be similar, and may be of any dimension. For example, we have also studied a four-dimensional system in which / 
and g are Henon maps with different parameters. For ease of presentation we restrict discussion to one-dimensional / 
and g. Systems such as Equation ([!]) are known in the mathematical literature as skew products or extensions; ours 
is constructed such that for c = 1, the x and y dynamics are identical and synchronized, whereas for c = 0, the x and 
y dynamics are completely independent. We investigate this system as c is decreased from 1 to 0. 

The simplest case is when f = g, for which the synchronization manifold M. (i.e. the line x — y) is invariant and 
attracting at c = 1. The bubbling bifurcation occurs when an orbit in M. loses transverse stability, typically via a 
period-doubling (pitchfork) bifurcation. This leads to the creation of new orbits outside of M. As the coupling is 
further reduced, more and more periodic orbits embedded in M. lose their transverse stability in a similar fashion 
, leading to the creation of additional orbits. As this process proceeds, the external UPOs simultaneously undergo 
period-doubling cascades to chaos, thus creating even more new orbits. We call the set of new orbits created in this 
fashion the emergent set. 

In the more general case / ^ g, x = y is by construction invariant and attracting for c = 1. Upon decreasing 
c, x = y is no longer invariant, and we observe that the UPOs migrate and spread out as shown in Figure |l| for 
coupled quadratic maps. As the coupling is decreased, we first observe transverse Cantor-like structure, followed by a 
"fattening" of the striations as the Lyapunov dimension of the attractor increases to 2.0 fig] . We have also observed 
similar UPO migration in the invertible case of coupled Henon maps. It is remarkable that this UPO migration 
appears to occur well before any orbit loses its transverse stability. In fact, we observe a large range of c over which, 
despite the apparent structure (ultimately two-dimensional), the periodic orbits migrate but are still transversely 
stable and one-to-one in the following sense: if the driver dynamics is fixed onto any one of its intrinsic period p 
orbits, then the limiting y dynamics is an attracting orbit of the same period. 

Let U be the set of unstable periodic orbits on the line x = y when c = 1. The number of orbits in U is determined 
by the driver and remains constant for all c because of the unidirectional coupling. For / = g, they remain fixed in 
place, but for f ^ g, they migrate as described above. As c is decreased from 1, the orbits' stability properties evolve, 
but they remain transversely attracting until a bubbling-type bifurcation is encountered. (We extend the concept 
of bubbling to the asymmetric case / ^ g by defining it as the point where the first orbit in U loses stability Jig].) 
As the coupling is further decreased, more and more orbits bifurcate and create orbits outside of U, and the above 
mechanism for the creation of the emergent set applies. Because of their migration, however, the orbits of U become 
intermingled among those of the emergent set. 

We wish to view system (Q) as an infinite collection of subsystems defined as follows. First, enumerate the periodic 
orbits of / (the driver), assigning each an index i = 1,2,.... Then subsystem Si is given by Equation (Q), but with the 
driver dynamics / locked on orbit i. The bifurcations described above correspond to bifurcations in the y components 
of these subsystems. Indeed, each subsystem Si exhibits a complete bifurcation structure in y as c is varied from 1 to 
0. 

To quantify the discussion, let N xy {p) denote the number of period p orbits of system (||). Also, let Nf(p) denote 
the number of period p solutions of f p {x) — x — alone; N g (p) is defined analogously. These quantities contain 
contributions from all periodic orbits of period q, where q is an integer factor of p. For c = 1, the system exhibits 
identical synchrony and N xy (p) = Nf(p). In contrast, when c = 0, the system is fully decoupled into independent 
systems, and N xy (p) = N f{p)N g {p). (Note that N xy (p) may achieve its maximum at c = or at intermediate values 
of c, depending on the nature of g.) 

Our goal is to elucidate how this change in the unstable periodic orbit structure proceeds as c is varied. To this end, 
we consider the topological entropy h [|l7| ; for large p, the number of periodic orbits of period p in a chaotic set increases 
exponentially withp as N(p) ~ e hp . Thus, the topological entropy of the coupled system is h xy = linip^oo lnN xy (p)/p, 
and similarly, the topological entropy of the driver is hf = linip^oo InN f(p) /p. 

Let N e (p) be the number of periodic orbits of period p that are not in U. These orbits reside in the emergent set 
and are created by the bifurcations described above. Thus N e (p) = Ylili Ni(p), where the summation is over the 
number rif, of subsystems that have bifurcated. N{(p) is the number of periodic orbits of period p not in U that are 
associated with a particular subsystem Si in the summation. The topological entropy of the emergent set is therefore 
h e = lim p ^ocNe(p)/p flf]- 

We define the state of topological coherence for system (|l|) as the condition h xy = hf. In this state, the topological 
entropy of the system is determined by the driver. In order for topological coherence to be destroyed, the topological 
entropy of the full system h xy must exceed h / . This occurs when the emergent set becomes sufficiently complex. In 
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the present case, h xy = linip^oo (Nf(p) + N e (p)) = max(hf,h e ). There is therefore a critical value a of coupling 
where h e first exceeds hf and topological coherence is lost. We call this the decoherence transition. For the special 
case / = g, we find that this typically occurs between the bubbling Cb u and the blowout Cbo bifurcations. 

We now address the question of how the decoherence transition may be measured from experimental data. First, we 
observe that in the symmetric noise- free case (/ = g), trajectories collapse onto M. and remain there until the blowout 
bifurcation. Thus, estimates of the decoherence transition based on measured data will not reflect the contribution 
of the emergent set. However, this case is exceptional. In the more general asymmetric case (/ ^ g), the orbits 
of U migrate and become intermingled with those of the emergent set. Because of this, typical trajectories do not 
necessarily remain near U; instead, the observed attractor incorporates parts of the emergent set. How much of the 
emergent set is incorporated depends on the degree of asymmetry and the coupling. By using trajectory data, an 
effective decoherence transition can be measured which indicates how much the emergent set actually influences the 
observed dynamics. It is precisely this effective transition that is most relevant to the observed dynamics of the 
system and hence is most relevant to experimental situations. Below we describe an efficient method for estimating 
the effective decoherence transition that is based on actual trajectory information. 

We use the methods of Ref. These authors define an average n-step stretching rate as follows. Let denote 
the square root of the largest eigenvalue of [J™(xi)] T J n (xj) for some initial condition x^, where J is the Jacobian of 
the system in Equation (TP). Then form the following average quantity over m initial conditions chosen with respect 

to the natural measure: ln{U) n , m = In [yZx W /™) For hyperbolic systems with one stretching direction, it 
can be shown ]l7fl that h\ — lim n _Kx) limm-too ln{ld) nim is the topological entropy of the system JTsj ] . Numerically 
hi can be obtained by measuring the scaling of InQ^ILx ^h^) with n, for a sufficiently large m. 

Since Equation ([!]) can exhibit two expanding directions for certain parameter values, we also measure the two 

dimensional average stretching rate: in(2d)„ )TO = In fX)^Li(^i^2)j™V TO J / n J where {Xi^)^ represents the square 

root of the product of the two largest eigenvalues of [J n (xj)] T J™(xj). The topological entropy when there are two 
expanding directions is given by hi = lim ra ^oo linim—xx, ln(2d) n , m . 

These quantities enable us to calculate the topological entropy of the full system as the coupling parameter varies 
and traverses regions with one and two stretching directions: h xy = max(hi, hi)- For the examples considered below, 
we have h\ = hf for the entire coupling range of interest. Thus, the effective decoherence transition occurs when hi 
first exceeds h\ = hf. (In higher dimensional systems, these methods are practical if the number of unstable directions 
is low; otherwise, it may be difficult to accurately calculate these quantities from limited data 20 1.) 

We apply these methods to a system of coupled quadratic maps. We take f(x) = 1.7 — x 2 , g(y) — a g — y 2 , and 
consider the cases a g — 2.0, 1.7, and 1.5. (We have also used two Henon maps coupled as in Eq. ffl, and the results 
are qualitatively the same.) Figure || shows max(hi, hi) for these cases. In all cases, h\ is equal to the topological 
entropy of the driver dynamics hf (for a g = 2.0, hi > hf only for c < 0.1, not shown). The effective decoherence 
transition occurs when hi exceeds h±, as indicated by arrows. 

Finally, we illustrate the influence of the emergent set for the symmetric case a/ = a g with noise. As c is decreased 
from bubbling, the dynamics will make occasional excursions from M.. We expect progressively longer transient times 
outside of M. due to the increasing complexity of the emergent set. To illustrate, we consider the symmetric case 
a f = a 9 = 1-7 and plot in Figure || the average duration of a burst versus coupling value. We define a burst as an 
excursion of at least 10 iterations beyond a small distance i5 = 0.05 from the line x = y\ the end is when the orbit 
falls back within this distance. We measure iterations per burst over trajectories of 10 6 iterations, then average over 
100 such realizations. For each noise level (uniformly distributed noise within a given amplitude), there is a transition 
range of c, depending on S and the magnitude of the noise, above which bursts are not observed. For very small noise, 
this transition is close to the blowout bifurcation; for larger noise, the transition shifts to higher values of coupling. 
Note that for c values below their respective transitions, the various curves asymptote to a common curve, suggesting 
that the dynamics during the bursts is consistently influenced by the emergent set. As expected, the average duration 
of bursts increases with decreasing c. 

Finally, we note that the mechanism for the creation of emergent sets outlined here leads one naturally to expect 
unstable dimension variability ||l| as a typical feature of emergent sets, and of coupled systems in general. 

In conclusion, we emphasize that the emergent set framework developed here is quite general and applies to coupled 
systems of non-identical elements where previously studied bifurcation structures may be problematic or inappropriate. 
Furthermore, the effective decoherence transition can be estimated in such systems from experimental data. 

We thank C. Grebogi, E. Ott, B. Hunt, and J. A. Yorke. This work was supported by NSF (IBN 9727739, E.B. 
and P.S.) and NIH (7K02MH01493, S.J.S.; 2R01MH50006, S.J.S., P.S., and B.J.G.). Figure @ was generated using 
DYNAMICS §1. 
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FIG. 1. Magnifications of the attractor for Equation with f(x) = 1.7 — x 2 , g(y) = 1.5 — y 2 , and c = 0.45 (inset, c = 0.9). 
Superimposed on both views are periodic orbits (+) of periods 1 to 20. At c = 1, these orbits lie on the diagonal x — y, but 
as c is decreased, they migrate as shown here. The Lyapunov dimension is 2.0 (inset, 1.25), and the orbits in both cases are 
transversely stable and one-to-one in the sense described in the text. 



FIG. 2. Topological entropy estimate max{h\, hi) for the system in Equation ([j]) with f(x) = 1.7 — x 2 , g(y) — a g — y 2 , and 
cases a g — 2.0, 1.7, and 1.5. hi in all cases equals hf, the topological entropy of the driver, for the range of coupling shown. 
The effective decoherence transition occurs when hi first exceeds h\ (arrows). (For the symmetric case a g = a,f = 1.7, bubbling 
occurs at ct u — 0.442 and blowout at c bo ~ 0.352.) 



FIG. 3. The duration of bursts away from the synchronization manifold for Equation (|l|) with f(x) — 1.7 — x 2 and 
g(y) — 1.7 — y 2 . The threshold value of c for observing such bursts increases with larger noise. The curves asymptote to a 
common curve for coupling less than this threshold, indicating that the trajectory during the bursts is strongly influenced by 
the emergent set. 
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Estimation of Effective Decoherence Transition 
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